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We investigate the entanglement features of the quantum states employed in quantum algorithms. 
In particular, we analyse the multipartite entanglement properties in the Deutsch-Jozsa, Grover and 
Simon algorithms. Our results show that for these algorithms most instances involve multipartite 
entanglement. 
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Entanglement is a major resource in quantum infor- 
mation processing. However, its role in achieving the 
quantum computational speedup in the currently known 
quantum algorithms is not yet completely clear and has 
been a highly debated question since the advent of quan- 
tum computation. In particular, the role of genuine 
multipartite entanglement has not been elucidated thor- 
oughly. It was shown that in Shor's algorithm multi- 
partite entanglement is needed to achieve exponential 
computational speedup with quantum resources |il]. In 
this work we analyse other known quantum algorithms, 
namely Deutsch-Jozsa, Grover and Simon's, and assess 
the multipartite entanglement properties of the pure 
quantum states employed. 

Let us first focus on the Deutsch-Jozsa [21] and Grover 
[3| algorithms, following the formulation given in [J]. The 
n-qubit states that occur in these algorithms are of the 
form 
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where \x) represent the computational basis states of n 
qubits and f{x) is the {0, 1}" — >■ {0, 1} Boolean function 
that needs to be evaluated. Notice that (— l)-''^^' = ±1 
is just a real phase factor. 

The above states are achieved by starting from the 
equally weighted superposition of all possible 2" states 
in the computational basis, namely the state 
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and then applying a unitary transformation Uf which 
involves an additional qubit (a target qubit) and acts as 
follows 



Uf\x)\y)^\x)\f{x)®y). 



(3) 



In the above expression \y) is the state of the target qubit 
and ® denotes addition modulo 2. The target qubit is 
initially prepared in the state |— ) = (|0) — |l))/\/2- In 



this way the target state is left unchanged by the action 
of Uf and the state of the n-qubit register takes the form 
([1]). We will now discuss the entanglement properties of 
the multi-qubit "real equally weighted states" lipf). 

In the Deutsch-Jozsa algorithm the function / to be 
evaluated is promised to be either constant or balanced 
(balanced means that it takes output for half of the 
inputs and output 1 for the others). The state ([T]), where 
a minus sign is present in front of half of the computa- 
tional basis states and will therefore be called balanced 
state, can be either separable or entangled. The state ([T]) 
is separable if and only if it can be expressed in the form 
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where a is an n-bit string and a-x = aixi(Ba2X2(B...anXn- 
Therefore, apart from global phase factors, the number 
of dinstinct separable states of the form ([T]) is 2". If the 
function / is constant the state {ipf) is clearly separable, 
and apart from a global phase factor, it is the state \ipo) 
(in this case the string a in the form ^ has all bits 
0). Notice that the above state in Eq. ^ corresponds 
exactly to the state computed in the Bernstein- Vazirani 
problem [5| , where the function is promised to be of the 
form fa{x) = a-x and the question is to find a. Therefore, 
for the Bernstein- Vazirani problem, which can be solved 
as the Deutsch-Jozsa problem restricted to the particular 
class of functions of the form /a, the qubits involved are 
never entangled. This has already been observed in [6|]. 
In general, if the function / is balanced the state ([T]) 
can then be either separable or entangled. We will now 
evaluate the number of separable versus entangled states. 
The number of possible balanced functions is given by 
Nbai = i?(2"', 2"'"^), where B denotes the binomial co- 
efficient. The number of balanced functions correspond- 
ing to separable states, which are then necessarily of the 
form (j4|), is Nbai,sep = 2(2" — 1). This value is obtained 
by noticing that two different balanced functions corre- 
spond to each a (they would correspond to states |i/)s) and 
— I'f/'s), which differ just in a global phase factor), while 
the case of vanishing a does not correspond to a balanced 
function and therefore must be subtracted. The fraction 
of balanced functions that involve separable states with 



respect to the total is therefore given by Nbai,8ep/Nbai = 
2(2" - l)(2"-M)2/2"!. Using the Stirhng approximation, 
i.e. lim2:^oo 2^! ~ \/2Trxx^ e"'' , we arrive at the asymp- 
totic hmit lim^^oo Nbai,sep/Nbai « %/2^(2" - l)2t/22". 
Thus, for large numbers of qubits the fraction of sep- 
arable balanced states becomes exponentially small, i.e. 
for most balanced functions the Deutsch-Jozsa algorithm 
employs entangled states. 

In the following we will denote by Sq the set of pure q- 
separable states, i.e. states that can be written as tensor 
products of pure states of q subsystems [7|]. So far we 
have evaluated the fraction of fully separable states (i.e. 
G Sn)- We can also evaluate the fraction of states that 
are not fully separable but contain entanglement only 
between two qubits (i.e. G 5'„_i\5„). Actually, this 
kind of states are of the form \^psep,n-2) 'S^ \'4'ent,2)- There 
are 2(2"~^ — 1) * 8 dinstinct balanced functions which 
correspond to states of this form: the factor 2(2"^^ — 1) is 
the number of balanced functions corresponding to fully 
separable states of n — 2 qubits, while 8 is the number 
of entangled real equally weighted states of two qubits 
(four states with just a minus sign in the superposition 
and four states with three minus signs). States of this 
type can occur for B{n, 2) different partitions of the n 
qubits, therefore the total number of balanced functions 
corresponding to states in iS'ji_i\S'„ is 2(2"~^ — 1) * 8 * 
B{n,2). By taking the limit of large n we can see that 
the fraction of these functions over the total number Nbai 
is still exponentially small. 

Rather than evaluating all possible classes of states 
we will answer the question whether most states in the 
Deutsch-Jozsa algorithm (in the limit n — )> oo) are gen- 
uinely multipartite entangled. To this end, we count the 
number of biseparable balanced states, i.e. states in S2 
corresponding to balanced functions that are of the form 



llpbisep) = \i>k) (8) \i>n-k) 



(5) 



and their fraction within all balanced states. Here, it 
does not matter whether the constituting states are en- 
tangled or separable, because any pure state that cannot 
be written in a biseparable form is fully entangled. Be- 
fore counting the number of biseparable balanced states, 
we will establish when a biseparable state is balanced. 

Lemma: For a pure real equally weighted state of n 
qubits which is g-separable, i.e. 

\ipq-sep) = |^fci)'8)|V'fe2)<8'--#fcp), with y^^kj^n; q> 1 

(6) 
where each |'0fci) is real equally weighted, the state 
IV'g-sep) is balanced if and only if at least one \ipki) is 
balanced. 

Proof: Denote by N^,_^ the number of plus (minus) 
signs of a real equally weighted state with k qubits. 
"=>": For a balanced n-qubit state we have N^' = N_'. 



For each partition with ki qubits in one subsystem, and 

(k) 

n — ki in the other subsystem, we thus have 7V_j_ • 
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solution of this equation is either N_ 



(n-ki) 
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N_ ' or 



ik'] (k) 

N^ — N_ , i.e. at least one of the two subsystems is 

in a balanced state. This argument holds for all possible 

partitions. 

"<^": Assume without loss of generality that \ipki) is 



balanced, i.e. N 
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The number of plus 
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signs in the n-qubit state is A^, 

N_ • N_ , and the number of total minus signs 



is N, 

Aki 
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^(*:i) . j\j-{n-ki) 



^(ki) . n["-'''\ Due to 



Ai'^^' we find Afi") = A^i"\ 



Al"-'' = Nr" we find Af^"' = A^i"^ D 

To count all biseparable balanced states we first fix fc, 
and also fix the partition in Eq. ^. The number of real 
equally weighted states where at least one of the two sub- 
systems is balanced is Nbisep{k) = 5(2*-', 2'°"^) • 2^" -|- 
B(2"-'',2"-''-i) • 2"^" - B(2'=,2'='^i) • B(2""'',2""*="i). 
This expression is derived by counting the number of bal- 
anced functions in the left term times the number of all 
functions in the right, plus vice versa, minus the terms 
where both parts are balanced and we have already in- 
cluded it before. Next, we have to sum over all possi- 
ble bipartitions for fixed k (which leads to the binomial 
B{n^ k) as factor) and then have to sum over all k. The 
only partition where this argument does not hold is the 
case oik = n/2 for even n. Here the factor 1/2 is needed 
to ensure that the partitions are not counted twice. Thus, 
the number of biseparable balanced states is given by 



L("-i)/2j 
KLp = E B{n,k)Nusep{k) 

k=l 
+ 2-^^"' "/2)^b<sep(n/2)(5„/2_L„/2j,0 



(7) 



We now want to find livan^oo N^^ /Ntai- To this end, 
it is convenient to write the above expression in the fol- 
lowing form 



Ni^ip = 5]S(n,fc)(B(2^2^-l).22"-^ (8) 

-5(2'=, 2'=-!) • B(2"-^2"-'=-l)) . 
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For n — >■ 00 the largest term in the above equation is the 
one with fc = 1, and therefore we find that 



hm A^,^4< hm 2(n-l)n.22"- . (9) 

n^-oo ^ n— >oo 



Using this upper bound, we arrive at 



atDJ 
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(10) 



Thus, in the limit of n — > oo the number of biseparable 
states among the balanced ones goes to zero, and we 
conclude that for large n the Deutsch-Jozsa algorithm 
typically employs genuine multipartite entanglement. 

We will now discuss the case of Grover's algorithm. 
The state ^ is achieved after the first application of the 
oracle. In this case the function / has output 1 for en- 
tries X that correspond to solutions of the search problem 
and output for values of x that are not solutions. Let 
us denote with M the number of solutions, which typ- 
ically is much smaller than the total number of entries 
2". The number of possible states of the form \4>f) is 
7Vm = B(2",A/). 

In the case of a single solution M — 1 the state ([1]) 
corresponds to an equally weighted superposition of all 
possible computational basis states with the same rela- 
tive phases, except for a single term, which has relative 
phase -1 with respect to the others. Such a state is fully 
entangled, because a biseparable state would contain at 
least two minus signs. Therefore, it is genuine multi- 
partite entangled. The number of possible states with 
Af = 1 is clearly 2", because the -1 phase can be in front 
of any of the 2" computational basis states. 

In general, we remember that a necessary condition for 
the state ([T]) to be fully separable (i.e. of the form ^), is 
that the output of / is constant or balanced. Therefore, 
for non-balanced or constant functions, namely whenever 
M ^ 0, 2"~^, the state ([T]) is always entangled. Moreover, 
in the case of odd M it is never possible to write the state 
as a tensor product because this would lead to an even 
number of -1 relative phases in the state IV'/)- Therefore, 
whenever the number of solutions M is odd the state \ipf) 
is fully entangled. 

Let us now consider the case M = 2. Here the state 
([T]) can be either biseparable or entangled. Actually, the 
only possibility to write it as tensor product would be in 
the form 
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(|O)+|l))0|V'e«t,«-l) , (11) 



where \ipent,n-i) represents a fully entangled state of n— 1 
qubits, of the form ([1]) and with just one relative phase 
-1 in the superposition, and the partition is arbitrary. If 
the state cannot be written in this form in any partition, 
then it is fully entangled. As mentioned above, the to- 
tal number of possible states with M = 2 is 5(2", 2), 
while the number of biseparable states among these is 
n2"~^. Therefore, the fraction of biseparable states for 



Af = 2 is exponentially small in the asymptotic limit. In 
other words, for large n and M = 2 the Grover algorithm 
typically employs genuine multipartite entanglement. 

Let us now consider the case Af = 4. The state 1-0/) in 
this case is biseparable, triseparable or fully entangled. 
Actually, it could be factorised in the triseparable form 
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where \ipent,n-2) represents a fully entangled state of ti— 2 
qubits, of the form ^ and with just one relative phase -1 
in the superposition, and the partition is again arbitrary. 
The number of possible triseparable states of the above 
form is 2"~^f3(n, 2). The only other possibility to factor 
it is in the biseparable form ([Tl]) . where \'ipent,n-i) repre- 
sents a fully entangled state of n — 1 qubits, of the form 
([1]) and with two components with -1 relative phases in 
the superposition. If the state cannot be written in the 
above biseparable form or in the triseparable form ([12]) 
for any partition, then it is fully entangled. 

The total number of biseparable states for M — 4 
(written as tensor product either in the triseparable form 
dni) or in the form (HH)) is iV^^^^p = nB(2"-\2). The 
fraction of biseparable states is then given by 



lim 
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nJ5(2"-\2) 
B(2",4) 



, (13) 



namely the states employed are typically genuine multi- 
partite entangled. 

The above argument can be generalised to values of M 
which are powers of 2: for Af = 2^ , the initial state |V'/) 
can be biseparable, triseparable, ... {k -f l)-separable or 
fully entangled. When the state is fc -I- 1-separable it has 
to be of the form 



IV'fc+l-sep) 



2fe/2 



(|0) + |l))^'=®|V'e„t,„-fc) 



(14) 



where \'4>ent,n.-k) is a fully entangled state of n — A: 
qubits with a single term in the superposition having 
relative phase -1. The reason for this structure is as fol- 
lows: assume that the k parties would contain entan- 
glement; then the number M of solutions, according to 
the Lemma, would be at least Af > 2"/^, while we re- 
quire Af <C 2". The number of states of the form (|14p is 
2"~'=B(n,A:). 

When the state is fc-separable it has to be of the form 



1^, 



k-sep! — 2(fe-l)/2 



m + \l)T''-^®\^ent,n-k+l) . (15) 



where now \ij}ent,n-k+i) is a fully entangled state of 
n ~ k + \ qubits with two components having -1 rela- 
tive phase. In general, for M = 2^ , the state 1-0/) can be 
j-separable, with j < fc -I- 1, or fully entangled. 



For a generic even M = 2'^{2p + 1), the state \ijjf) can 
be always biseparable or fully entangled. In particular, 
the state can then be j-separable, with all values of j 
ranging from 2 to g + 1, or fully entangled. The number 
of biseparable states for even M is in general given by 
^bisepi^) = nB{2"-\M/2). Therefore, the fraction of 
biseparable states is given by 
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N^^sepiM) 
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M 



= lim ^— ' / ' ~ , (16 



where we consider q and p finite and fixed. We can then 
conclude that for any M ^ N (as mentioned above, in 
our analysis it is sufficient that M < 2"/^) the states are 
typically multipartite entangled. 

We will now consider Simon's algorithm [Sj. In this 
case the function to be evaluated is promised to be a 
{0, 1}" -> {0, 1}" periodic 2 — >• 1 function and the task is 
to evaluate the period r (r is an n-bit string, with r ^ 0) 
with the smallest number of evaluations. In other words, 
in Simon's case f{x) = f{y) if and only ii x = y®r. The 
first register, as in the previous cases, is composed of n 
qubits prepared in state IV'o); while the target register is 
now composed of another set of n qubits, which are all 
prepared in state |0). 

After the function evaluation step ([3]), the global state 
of the two registers takes the form 



i41„j^^==y:(ix,) + |x, 
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r))\fix,)), (17) 



where Xi is the set of 2"^^ input values leading to dif- 
ferent outputs. Notice that the above state is always 
entangled between the two registers. The next step is 
a measurement on the second register, after which the 
state of the first register takes the simple form 



|Vs,™o„> = -^(|x) + |5 + r)) , (18) 

where x is now a random value among the Xi^s. We will 
now study the entanglement properties of the states of 
the form (IT5|) . Notice that, by keeping r fixed and chang- 
ing X in the above form, the entanglement properties of 
the state do not change because by applying just local 
operations (of the type a^) we can reach all the states of 
the above form (with different x and same r). We can 
then for simplicity analyse the properties of these states 
by choosing the particular value x = 0, namely 



vanishing values of r, is 2" — 1. If r contains a single 
value 1, i.e. wt{r) = 1 {wt{r) is the weight of the binary 
string r, namely the number of I's that it contains) then 
the state |V'Simon,o) is fully separable (g S^), because it 
isoftheforml/\/2(|0)-F|l))|0)...|0). There are n possible 
states of this type, corresponding to the different strings 
of length n with weight 1. 

If wt{r) = 2, then the state (fT9|) contains two-qubits 
entanglement (g 5'„_i). The number of states of this 
form is B{n,2). This simple counting argument can be 
generalised to the case of an arbitrary weight for r: if 
wt{r) = k the state P^ belongs to Sn-k+i and there 
are B{n^ k) of them. The case of maximum weight is 
wt{k) = n and then the state is n-partite entangled, of 
the GHZ form. Notice that the most populated class of 
states corresponds to values of the period r with wt{r) — 
[n/2j , i.e. typically half of the qubits are entangled with 
each other. 



In summary, we have elucidated the role of multi- 
partite entanglement in the Deutsch-Jozsa, the Grover 
and the Simon algorithms. For the Deutsch-Jozsa al- 
gorithm, multipartite entanglement within the first reg- 
ister is needed to accommodate all possible (balanced) 
functions. We have shown that the fraction of balanced 
functions corresponding to biseparable states decreases 
exponentially with the number of qubits. In this sense 
most balanced functions involve genuine multipartite en- 
tanglement. In the Grover algorithm, we have shown 
that the entanglement properties of the initial state of 
the first register depend on the number of solutions and 
we have demonstrated that such a state is also typically 
multipartite entangled, when a small number of items 
is searched for. In Simon's algorithm the situation is 
slightly different: in general the first register and the tar- 
get register are entangled. Moreover, the kind of entan- 
glement within the first register after the measurement 
of the target register depends on the weight of the period 
r: for increasing weight it involves an increasing number 
of qubits. For large n, the states employed typically in- 
volve entanglement among half of the constituent qubits. 
We have thus shown that multipartite entanglement is an 
essential property in the considered quantum algorithms. 
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|V'S^mon.o) = -/|(|0> + |r)) 



(19) 



Notice that the number of possible different states of 
the above form, corresponding to all the possible non 
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